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Abstract 

The most advanced implementation of adaptive constraint processing with Constraint 
Handling Rules (CHR) allows the application of intelligent search strategies to solve Con- 
straint Satisfaction Problems (CSP). This presentation compares an improved version of 
conflict-directed backjumping and two variants of dynamic backtracking with respect to 
chronological backtracking on some of the AIM instances which are a benchmark set of 
random 3-SAT problems. A CHR implementation of a Boolean constraint solver com- 
bined with these different search strategies in Java is thus being compared with a CHR 
implementation of the same Boolean constraint solver combined with chronological back- 
tracking in SICStus Prolog. This comparison shows that the addition of "intelligence" to 
the search process may reduce the number of search steps dramatically. Furthermore, the 
runtime of their Java implementations is in most cases faster than the implementations of 
chronological backtracking. More specifically, conflict-directed backjumping is even faster 
than the SICStus Prolog implementation of chronological backtracking, although our Java 
implementation of CHR lacks the optimisations made in the SICStus Prolog system. 

KEYWORDS: dynamic backtracking, conflict-directed backjumping, rule-based constraint 
handling, intelligent search, SAT problems 



1 Introduction 

Constraint Handling Rules (CHR) are multiheaded, guarded rules used to propagate 
new or simplify given constraints (Fruhwirth 1995 Fruhwirth 1998). For example, 
the CHR 

leq(X.Y), leq(Y,Z) ==> leq(X.Z) . 

reflects the transitivity of the binary relation leq. Thus, for any two constraints 
leq(A,B) and leq(B,C) an additional constraint leq(A,C) is derived - implicitly 
given knowledge is made explicit. Another CHR 

leq(X,Y), leq(Y.X) <=> X=Y. 

reflects the symmetry of the binary relation leq. Thus, any two constraints 
leq(A,B) and leq(B,A) are replaced by the syntactical equation A=B. 
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A detailed formal description of the syntax, the declarative and operational se- 
mantics of CHR is omitted in this paper because these topics are addressed in depth 
in the literature, e.g. in IjFriihwirth 1 998). 

There are several CHR implementations, e.g. in 
ECLiPSe (IFrvihwirth and Brisset 1995fc . in SICStus Pro- 

log (jHolzbaur and Friihwirth 2000jl or even in Java IjSchmauss 19991 IWolf 2001ajl . 
All but the last IjWolf 2001a|) only support constraint deletions implicitly through 
chronological backtracking. Arbitrary sequences of constraint additions and 
deletions, which are necessary for intelligent search strategies like conflict-directed 
backjumping l|Prosser 19931 IPro sscr 199E>J) or dynamic backtracking (Bak er 19941 
|Ginsberg 1993| I.Tussien et al. 2000|) . are not supported. Furthermore, if there is an 
inconsistency, the "classical" CHR implementations offer users no help in finding 
out what causes this inconsistency. 

This paper reviews the first implementation of "adaptive" CHR, cf. ( Wol f 19991 
IWolf 2001aJ) . In this context, "adaptive" means that constraint additions and dele- 
tions in arbitrary order are supported, i.e. after each change of the considered 
constraints, the derivations based on CHR are adapted accordingly. Thus, deletion 
of the constraint leq(A,B) or leq(B,C) causes the derived constraint leq(A,C) to 
be deleted, too. However, this implementation is in Java, which does not support 
backtracking like Prolog systems; thus depth-first search is not intrinsic, enabling 
different search strategies to be realized directly and not on top of the underlying 
chronological backtracking mechanism. Additionally, this implementation returns 
an explanation for any occurring inconsistency, thus supporting explanation-based 
constraint programming ( Juss ien 200 l|l . This allows not only user guidance, e.g. 
during debugging of incorrect constraint models or in interactive constraint solving, 
but also automatic constraint relaxation as well as dynamic problem handling in 
reactive systems. Moreover, explanations can be used to build new "explanation- 
directed" search algorithms. 

The aim of the paper is to show that this adaptive CHR implementation is very 
well suited for implementing not only depth-first search but also intelligent search 
algorithms like sophisticated conflict-directed backjumping and dynamic backtrack- 
ing, while maintaining consistency. The given implementations show that 

• constraint propagation ideally replaces the proposed constraint checks/tests 
in these intelligent search algorithms 

• justifications of (derived) constraints, especially of false, properly act as 
conflict sets in conflict-directed backjumping or as elimination explanations 
in dynamic backtracking 

• the possibility of arbitrary constraint deletions directly supports non- 
chronological backtracking 

• constraint handling (i.e. propagation) maintains (local) consistency, offering 
early detection and good avoidance of dead ends during the search 

To be more precise, the implementation of these algorithms is specialised for 
Boolean constraint problems where the variables have only two possible values: 
and 1. However, a generalisation for other (finite) domains is quite simple because 
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the interaction with the Boolean constraint solver written in CHR is opaque. We 
therefore assume that any other terminating constraint solver realized with CHR 
will work as well. We cite the soundness and completeness of CHR (Fruh wTrth 1998|) 
as well as the correctness and termination of the adaptation of CHR deriva- 
tions I) Wolf 1999)) based on explanations as evidence for this claim. 

The paper is organised as follows. The next section briefly looks at the adaptive 
CHR system. Section [3] presents a CHR-based specification of a Boolean constraint 
solver to solve SAT(isfiability) problems formulated as prepositional logic formulas. 
The compilation process of these rules into an adaptive constraint solver is explained 
and the application programming interface (API) for this solver is described. Sec- 
tion 01 introduces the AIM instances, a benchmark set of random 3-SAT problems 
containing instances with exactly one solution and instances that are inconsistent. 
Section presents our implementations of different search strategies, from chrono- 
logical backtracking fSection l5.3|) to conflict-directed backjumping (Section El} to 
dynamic backtracking fSection l5.5(l . These implementations are built on top of the 
Boolean constraint solver presented in Section |2| These solvers are applied to all 
AIM instances with 50 variables, either to solve them or to detect their inconsis- 
tency. Section shows and compares the required backtracking/backjumping steps 
and their runtime. Section attempts to analyse the measured results. Section |H1 
concludes the paper. 



2 The Adaptive CHR System 




generates 



Initially, the adaptive CHR system consists of a 
runtime system and a compiler. They contain the 
data structures that are required to generate rule- 
based adaptive constraint solvers and to implement 
Java programs that apply these solvers to dynamic 
CSP. The definition of a rule-based constraint solver 
is quite simple: the CHR that define the solver for 
a specific domain are coded by the user in a so- 
called CHR handler. Here, a CHR handler consists 
of Java objects representing CHR which are com- 
piled to Java programs by the use of the compiler. 
Compiling and running a CHR handler generates a 
Java package containing Java code that implements 
the defined solver and its interface: the addition or 
deletion of user-defined constraints or syntactical equations, a consistency test and 
the explanation of inconsistencies. These methods allow dynamic constraint solving 
as well as explanation-based constraint programming ( Juss ien 200 If) in any applica- 
tion: 

• Constraints may be added and deleted in arbitrary order. 

• Constraint handling, i.e. propagation, is performed accordingly. 

• Whenever an inconsistency is detected, the explanation identifies a subset of 
constraints causing this inconsistency. 



Fig. 1. The architecture of 
the adaptive CHR system 
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A user application interacts with the CHR package provided by the user in the 
CHR handler and the runtime system. Figure shows the components and their 
interactions. 

During compilation for each handler, a constraint system class is generated re- 
taining the name of the handler. Furthermore, for each head constraint of a CHR, 
a method of this class is generated retaining the name and arity of the CHR to add 
user-defined constraints to the constraint store. In addition to these handler-specific 
methods each constraint system class has common methods to justify the assign- 
ment of an integer to a variable, i.e. to add a syntactical equation justified by an 
integer to the constraint store; to delete all constraints with a specific justification, 

1. e. in a set of integers; to test the consistency of the currently valid syntactical 
equations; and to get an explanation, i.e. a set of justifications (integers) that is 
responsible for a detected inconsistency: 

• void equal (Variable var, int i, IntegerSet set) 

• void delete (IntegerSet set) 

• boolean isConsistent () 

• IntegerSet getExplanationO 

The class Variable implements logical variables that may be bound to logical 
terms (objects of the class Term), which are either numbers, logical variables or 
function terms. To simplify matters, integer sets (objects of the class IntegerSet) 
and operations on it are represented by the use of the usual mathematical set 
notation, e.g. the set consisting of the integers 2, 3, and 5 is represented by {2, 3, 5} 
and the union of two sets A and B is represented by A U B. 

Example 1 Let a CHR handler called trans consist of the CHR 

leq(X,Y), leq(Y.X) <=> X=Y. 

specifying the user-defined constraint leq. Furthermore, let the constraints 
leq(0,A) and leq(B, 1) with empty justifications be already added to the constraint 
store of the constraint system cs, i.e. be an object of the class trans. Assuming 
that A and B are constraint variables (objects of the class Variable/, cs. equal (A, 

2, {1}) adds the equation A=2 justified by the set 1 {1} to the constraint store of 
cs. Thus, the value 2 is assigned to the variable A and the store contains A=2, 
leq(0,2), both justified by {1}, and leq(B,l) justified by the empty set. Then, 
the call cs . isConsistent () returns true. Further addition of the equation B=2 
with the justification {3} is realized by calling cs. equal (B, 2, {3}). The result- 
ing constraint store now contains A=2, leq(0 ,2) , both justified by {1}, and B=2 and 
leq(2,l), both justified by {3}. This triggers the CHR, which replaces leq(0,2) 
and leq(2,l) by the equationO=l, i.e. an inconsistency justified by {1,3}. Thus, the 
call cs . isConsistent () returns false and cs . getExplanationO returns {1,3}. 
The detected inconsistency is eliminated by calling cs . delete ({1}). Afterwards, 
the constraint store contains leq(A,2) with the empty justification, and B=2 and 
leq(2,l) , both justified by {3}. □ 

1 For a unitised handling, integral identifiers are coded in singleton integer sets. 
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The next section contains a more relevant, practical example, illustrating how 
constraint solvers, i.e. CHR handlers, are defined and integrated into an application. 

3 A Rule-based Boolean Constraint Solver 

The ECLiPSe and SICStus Prolog distributions of CHR, or even WebCHR at 
http://www.pms.informatik.uni-muenchen.de/~webchr/, come with a simple 
but important constraint solver for Boolean constraints. This solver is essential 
for problems that are formulated as SAT problems, i.e. satisfiability problems of 
propositional logic formulas. The provided Boolean constraint solver supports the 
usual unary and binary operations on propositional variables: negation, conjunc- 
tion, disjunction (non-exclusive and exclusive) as well as implication. If we confine 
ourselves - without any loss of expressiveness - to problems in conjunctive normal 
form, only negation and disjunction have to be supported by a Boolean CHR solver 
as constraints, i.e. the disjunctions, are implicitly conjunctively connected by the 
separating comma. For instance, the formula in conjunctive normal form 

(A V ->B V C) A (-u4 V B V D) 

is equivalent to 

neg(A.F), neg(B.E), or(A,E,X), or(X,C,l), or(F,B,Y), or(Y,D,l) 

if the semantics of the user-defined constraint neg(X,Y) is ->X — Y and the se- 
mantics of the user-defined constraint or(X,Y,Z) is X V Y — Z for any arguments 
X, Y, and Z that are either propositional variables, 0, or 1. Thus, the important 
class of SAT problems may be modelled as constraint problems and solved by 
using the CHR handler with the following rules: 



or(0,X,Y) <=> Y=X. 
or(X,0,Y) <=> Y=X. 
or(X,Y,0) <=> X=0,Y=0. 
or(l,X,Y) <=> Y=l. 
or(X,l,Y) <=> Y=l. 
or(X,X,Z) <=> X=Z. 
neg(O.X) <=> X=l. 
neg(X.O) <=> X=l. 
neg(l.X) <=> X=0. 
neg(X.l) <=> X=0. 
neg(X.X) <=> fail. 



or(X,Y,A) \ or(X,Y,B) <=> A=B. 
or(X,Y,A) \ or(Y,X,B) <=> A=B. 
neg(X.Y) \ neg(Y.Z) <=> X=Z. 
neg(X.Y) \ neg(Z,Y) <=> X=Z. 
neg(Y.X) \ neg(Y.Z) <=> X=Z. 



neg(X.Y) \ or(X,Y,Z) 

neg(Y.X) \ or(X,Y,Z) 

neg(X.Z) , or(X,Y,Z) 

neg(Z.X) , or(X,Y,Z) 

neg(Y,Z) , or(X,Y,Z) 

neg(Z.Y) , or(X,Y,Z) 



<=> Z=i. 

<=> Z=l. 

<=> X=0,Y=i,Z=i. 

<=> X=0,Y=i,Z=i. 

<=> X=1,Y=0,Z=1. 

<=> X=i,Y=0,Z=i. 



The transformation of the CHR in the Java system (Wolf 2001a) is straightfor- 
ward: 



Example 2 The coding of the first rule or(0,X,Y) <=> Y=X. in a CHR handler 
is quite simple: 

class boolHandler { 

(01) public static void main(String [] args) { 

(02) DJCHR djchr = new DJCHR( "bool" , new String [] {"or/3" , "neg/2"}) ; 

(03) Variable x = new Variable ( "X" ) ; 
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(04) 
(05) 



Variable y = new Variable ( "Y" ) ; 
Term zero = new Term(0); 



(06) 



Term[] remove, keep, guard, body; 



(07) 
(08) 
(09) 



remove = new Term[]new Term("or" ,new Term [] {zero, x,y}) ; 

body = new Term[]{DJCHR.eq(y,x)}; 

dj chr . addRule (remove , null , body , null) ; 



(10) 

(11) } 
} 



dj chr . compileAllQ ; 



First of all, a new handler object djchr is generated (line 2). It is called bool 
and supports the ternary user-defined constraint or and the binary user-defined 
constraint neg. Then, two variables X and Y as well as a constant are gener- 
ated (lines 3-5). Every rule is split up into four arrays of terms (line 6): the head 
constraints that are removed, the head constraints that are kept, the guard con- 
straints, and the body constraints according to \Holzba,ur and, Fruhwirth 2000]) . For 
the considered rule to be transformed, the keep and guard arrays must be empty. 
However, the remove array contains the constraint or(0,X,Y), which is generated 
accordingly (line 7). Furthermore, the body constraint Y=X is generated using of the 
built-in method eq (line 8). Then the rule is composed and added to the handler 
(line 9). Finally, all added rules are compiled by calling compileAllO (line 10). □ 

During the compilation process, a Boolean constraint solver class called bool is 
generated. The application interface generated for this solver comprises the methods 

• void or_3(Term[] args, IntegerSet set) 

• void neq_2 (Term [] args, IntegerSet set) 

to add the specified user-defined constraints or and neg. 

Additionally, a class boolVariable of attributed logical variables, which is a sub- 
class of the class Variable, is also generated. It has special attributes to store and 
access efficiently the Boolean constraints on these variables, cf. iHolzbau r 19901 
IWolf 2001b(l . Thus, the propositional formula in conjunctive normal form 



is modelled as a Boolean constraint problem by the Java code fragment 

bool cs = new boolO; 
boolVariable a = new Variable ( "A" ) ; 
boolVariable b = new VariableC'B") ; 
boolVariable c = new VariableC'C") ; 
boolVariable d = new Variable ( "D" ) ; 
boolVariable e = new Variable ( "E" ) ; 
boolVariable f = new Variable ( "F" ) ; 
boolVariable x = new Variable ( "X" ) ; 
boolVariable y = new Variable ( "Y" ) ; 
Term one = new Term(l) ; 



(A V -lB V C) A (-.A V B V D) 
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cs.neg_2(new Term [] {a, f } , 0); 
cs. neg_2 (new Term []{b,e}, 0); 
cs.or_3(new Term[] {a,e ,x} , 0); 
cs.or_3(new Term [] {x, c , one} , 0); 
cs.or_3(new Term [] {f ,b,y} , 0); 
cs.or_3(new Term [] {y,d, one} , 0); 

if it is assumed that the formula is always valid, which means that the justifications 
are the empty sets. The calls of the methods cs .neg_2 and cs . or_3 add the con- 
straints to the constraint store of the constraint system cs and eventually trigger 
some of the compiled rules. 

It should be noted that the presented Boolean CHR solver applies the unit clause 
rule IjDavis and Putnam 1960|l . Unit clauses are disjunctions of literals, i.e. propo- 
sitional variables or their negations, where all literals except one are 0. Here, unit 
clauses are represented by conjunctions of k constraints 

or(X ,Xi,Ri), or(Ri,X 2 ,R2), • • • , or (R fe ,X k , 1) , 

where for a fixed index j 6 {1, . . . , k} it holds Xj = for all indices i ^ j. 

These constraints trigger the rule or(X,0,Y) <=> X=Y several times deriving in 
this order Rfc = . . . = Rj = 1, and further Ri = ... = Rj_i if j > 1 holds. In any 
case, either the rule or (X,0,Y) <=> X=Yor or(0,X,Y) <=> X=Y is finally triggered, 
which results in Xj = 1 in cither case. 

Other instances of prepositional formulas in conjunctive normal form that are 
processable using the introduced Boolean constraint solver are the AIM instances 
presented in the next section. 

4 The AIM Instances 

The AIM instances are random 3-SAT problem instances in conjunctive normal 
form, named after their originators Kazuo Iwama, Eiji Miyano and Yuichi Asahiro. 
3-SAT problems are conjunctions of disjunctions of three literals, i.e. propositional 
variables or negations of them. The AIM instances are all generated with a par- 
ticular random 3-SAT instance generator Ijlwama et al. 1996jl . The particularity is 
that the generator generates yes-instances and no-instances independently for wide 
ranges. Thus its primary role is to provide the sort of instances that conventional 
random generation has difficulty generating. The generator runs in a randomised 
fashion, which means that the 3-SAT instances essentially differ from those gener- 
ated in a deterministic fashion or from those translated from other problems. As a 
result, the following set of considered AIM instances includes 

• no- instances with low clause/ variable ratios that are inconsistent 

• yes-instances with low and high clause/ variable ratios that have exactly one 
solution 

The instances are called aim-xxx-y_y-zzzz-j where 

• xxx shows the number of variables, one of 50, 100 and 200 
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• y_y shows the clause/ variable ratio y.y, including 1.6, 2.0 for no-instances 
and 1.6, 2.0, 3.4, and 6.0 for single-solution yes-instances 

• zzzz is cither "no" or "yesl" , the former denoting a no-instance and the latter 
a single-solution yes-instance 

• the last j means simply the j-th instance at that parameter 

For each parameter, four instances are included in the benchmark set. The 
whole benchmark set is available online at http : //www . satlib . or g] For example, 
aim-50-l_6-no-l through aim-50-l_6-no-4 are four no- instances with 50 vari- 
ables and a 1.6 clause/variable ratio. In all, there are 18 sets of instances with 50, 
100 and 200 variables. For the yes-instances, clause /variable ratios are taken from 
1.6, 2.0, 3.4, and 6.0; for the no-instances, they are taken from 1.6, and 2.0. 

To find the unique solutions of the yes-instances or to prove the inconsistency 
of the no-instances, there are several state-of-the-art SAT solvers. A collection of 
SAT solvers is also available at http : //www . satlib . org Most of these algorithms 
are (heuristic) local-search algorithms or can be traced back to the Davis-Putnam 
procedure IjDavis and Putnam 1 960). However, the presented Boolean constraint 
solver in the previous section, complemented by a search procedure that assigns 
the value or 1 to the propositional variables, can obviously be used to solve such 
SAT problems. Furthermore, SAT problems are often used to compare "intelligent" 
search procedures with chronological backtracking, cf. ( Ba ker 19941 [Ginsberg 19931 
|Lynce and Marques-Silva 2002| IProsser T 993 ) . The next section therefore considers 
several search procedures and their interaction with the generated Boolean con- 
straint solver. 

5 The Search Procedures 

Before describing the compared search procedures in detail, we look at some of the 
assumptions made and programming conventions used. 

5.1 Programming Conventions 

It is assumed that there are Boolean Constraint Satisfaction Problems (CSP), i.e. 
there are variables Vi, ... ,V n with Boolean domains {0, 1}. Additionally, there are 
two types of Boolean constraints over these variables, either negations -^X = Y or 
disjunctions XVY — Z, where X, Y or Z are either variables, or 1. The problem 
is either to detect that there is no assignment of values to the variables such that 
the constraints are satisfied, i.e. the problem is inconsistent, or to find such an 
assignment, i.e. a solution. The Boolean constraints are realized by user-defined 
constraints handled by the Boolean solver presented in Section 01 

The different search procedures to solve Boolean CSP are presented in pseudo- 
code strongly related to Java. The main difference compared to Java is that math- 
ematical set notation is used instead of some methods of the "abstract" class 
IntegerSet. - Actually, we used our implementation of sparse integer sets, which 
is described in (|Wolf 1999J) . However, this might be replaced by any other, even 
more efficient implementation. 
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It is assumed that there is a globally declared array of variables 2 var, such that 
var [i] represents the variable Vi for i = 1, . . . , n where n is the actual number of 
variables in the considered problem. Variables (of the class Variable) implement 
attributed logical variables: they may be bound to terms, e.g. integers, or unbound, 
i.e. free. Thus, there is the method 

• boolean isBoundO which returns true if and only if the variable is bound. 

If a variable is bound, the method Term value () is defined, which returns the 
term the variable is bound to. Furthermore, there is an integer field num holding 
either the next value to be assigned to this variable (see Sect ion I5"3"jl or an identifier 
justifying the current assignment (see Section l5.5fl . 

A variable also contains an array of integer sets with indices ranging over the 
Boolean domain from to 1. If defined, i.e. if different from null, this array contains 
for each value unique identifiers of the variables, i.e. their indices, bound to values 
that result in an inconsistency, which was detected with respect to the considered 
Boolean constraint problem. 

Example 3 Let this set for the value 1 of the variable Vyj be {3, 7, 8, 10} where 
3,7,8, and 10 are the indices of other labelled variables. Then the assignment Vrj = 1 
is inconsistent with the current assignments to the variables Vs, V-j, Vs and Vio with 
respect to the considered Boolean CSP. □ 

In conflict-directed backjumping, these sets are called conflict sets, and in dy- 
namic backtracking they are called elimination explanations. Thus, in these search 
procedures the array is declared as either 

• IntegerSet [] conf lictSet or 

• IntegerSet [] elimExpl, 

accordingly. Furthermore, it is assumed that the language supports variable lists, 
e.g. a "wrapper" VariableList of the Java class ArrayList that supports 

• access to the size of a list: int sizeO 

• addition of a variable at the end of a list: void add (Variable var) 

• access to a variable at a specific position in a list: Variable get (int i), 
where the index of the first variable in a list is zero 

• access to the last variable in a list: Variable getLast(int i) 

• removal of a variable at a specific position: Variable remove (int i) 

such that the indices of the variables that come after the removed variable 
are decremented by one 

The following data structures are also assumed to be globally declared and thus 
accessible to all methods: 

• A unique Boolean constraint system cs of the class bool, where the con- 
straints are stored and processed by use of the Boolean CHR solver (see 
Section 



2 Variables in the sense of constraint processing. 
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static boolean cssp(int n) { 

(01) int i = 1; 

(02) while (1 <= i && i <= n) { 

(03) int j = xxxLabel(i); 

(04) if (i == j) 

(05) i = xxxUnlabel(i); 

(06) else i = j; 

(07) } 

(08) if (i = 0) 

(09) return false; 

(10) if (i > n) 

(11) return true; 

} 

Fig. 2. The cssp function for solving constraint satisfaction search problems 

• Two variable lists unlabelledVars and labelledVars of the class 
VariableList, where the unlabclled and labelled variables are stored dur- 
ing dynamic backtracking (cf. Figures ITUI and ITT1) . 

• A unique integer cntr, which is initially and incremented by one after an 
assignment in dynamic backtracking (cf. Figure I1UI line 10) serving as its 
unique justification. 

In the sequel, the calls to the constraint system using the interface to the adaptive 
CHR system are underlined. This shows the simple and powerful use of our adaptive 
CHR system in sophisticated search procedures. 

5.2 The Constraint Satisfaction Search Process 

According to the style presented in {Pressor 1993), the constraint satisfac- 
tion search problem (cssp) method in Figure establishes the environment 
in which the different search methods are called. The cssp method takes 
the total number of variables to be labelled with values and returns true if 
a solution is found and false if the given Boolean CSP is inconsistent. The 
"generic" methods xxxLabel and xxxUnlabel are replaced in the sequel re- 
sulting in chronological backtracking (cbtLabel/cbtUnlabel), conflict-directed 
backjumping (cbjLabel/cbjUnlabel) and two variants of dynamic backtracking 
(dbtLabel/dbtUnlabel and f btLabel/f btUnlabel). In all these instances, the 
method xxxLabel attempts to find a consistent assignment to the i-th variable. 3 
For this, the method takes i as its argument. It returns this given integer if no 
such assignment is found. However, if a consistent assignment to the i-th variable is 
found, it returns i + 1 after binding this variable to a value that is consistent with 
the other i — 1 previously bound variables and with respect to the given Boolean 
constraint problem. If xxxLabel returns i, the method xxxUnlabel is called. When 

3 The i-th variable coincides with V{ in chronological backtracking and conflict-directed backjump- 
ing but not necessarily in dynamic backtracking, which may dynamically change the variable 
ordering. 
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i + 1 is returned with 1 < i + 1 < n, xxxLabel is called again, looking for an 
assignment to the (i + l)-th variable. Returning n + 1 causes cssp to return true 
because a consistent assignment for all variables is found. 

The corresponding instance of xxxUnlabel is called when no consistent assign- 
ment to the i-th variable is found (cf. lines 4-5 in Figure^) ■ It performs backtracking 
from the i-th variable to an h-th variable (h < i) if another value for the h-th vari- 
able might resolve the inconsistency detected at the z-th variable. It takes i as its 
argument. It either returns or the index of the next variable to be labelled. Zero is 
returned if the detected inconsistency is not resolvable, i.e. the given Boolean CSP 
is inconsistent causing cssp to return false ( Figure|21 lines 8-9). 

5.3 Chronological Backtracking 

Chronological backtracking (CBT) is a simple 
depth-first search (cf. Figure |3J) with a fixed tree 
structure, i.e. variable ordering. If the variables are 
not already bound by constraint processing (Fig- 
ure 0] lines 1-2), they are incrementally bound to 
the values or 1. First, the current variable V$ is la- 
belled with the value 0. Search continues with Vi+i 
if no inconsistency is detected (Figure^] lines 5-6). 
Otherwise, the value 1 is assigned to the variable 
Vi. Again, search continues with Vi+i if no inconsis- Fig. 3. The principle of 
tency is detected. Otherwise, a dead end is reached chronological backtracking 
and the search backtracks to the variable VJ_i (cf. Figures El and EJ. 

int cbtLabel (int i) { 

(01) if (var[i] .isBoundO) 

(02) return i+1; 

(03) while (var[i] .num <= 1) { 

(04) cs . equal (var [i] , var [i] .num++, {i}) ; 

(05) if ( cs . isConsistent () ) 

(06) return i+1; 

(07) else 

(08) cs.delete({i}) ; 

(09) } 

(10) return i; 
} 

Fig. 4. The labelling method of chronological backtracking 

A generalisation of this search process for arbitrary finite domains is quite 
simple: the field num in the variable must be replaced by the domain. During 
labelling, it must be iterated over the values in the current domain (Figure 0] lines 
3-9). The iterator for this loop must be reset during unlabelling (Figure line 2). 
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int cbt Unlabel (int i) { 

(01) cs.delete({i}) ; 

(02) var [i] .num=0; 

(03) return i-i; 
} 



Fig. 5. The unlabclling method of chronological backtracking 




5-4 Conflict- directed Backjumping 

Conflict-directed backjumping (CBJ) (Prossc ~ 1993|l is a guided depth-first search 
with a fixed tree structure, i.e. variable ordering, that "jumps back" to the most 
recent variable assignment that is in conflict with the current variable (cf. Figure|HJ). 
Originally, CBJ maintains a conflict set per variable. However, in our refinement it 
maintains a conflict set for each value of every variable. Initially, these conflict sets 
are not defined, i.e. null. 

If the unlabelled variable V, is not already bound 
by constraint processing (Figure [7| lines 1-2) the 
attempt is made to bind it either to the value or 
1. The current variable Vi is labelled with the first 
value that is possibly not in conflict with other al- 
ready labelled variables (Figure [7| line 4-5). Thus, 
the index of the variable is chosen as the justi- 
fication of this assignment because it simply al- 
Fig. 6. The principle of \ ows arL y subsequent deletion of it and all its conse- 
conflict-directed backjump- quences computed by the underlying Boolean con- 
m S straint solver (cf. Figure [7\ line 10 and Figure [SJ 

line 6). 

The search continues with V^+i if no inconsistency is detected (Figure lines 
6-7). Otherwise, the indices of the already labelled variables that are responsible for 
the detected inconsistency form the conflict set of the attempted value (Figure [7| 
lines 8-11), the assignment is deleted (Figured line 10) and the next value for Vi 
is attempted (Figure lines 3-13). If all assignments lead to an inconsistency, a 
dead end is reached, i.e. i is returned (Figure[7| line 14), which triggers unlabelling. 
If the conflict sets of all values of the considered variable Vi are empty, the deletion 
of all variable assignments will not resolve the detected inconsistency. Thus, is 
returned, indicating the inconsistency of the given Boolean CSP (Figure 03 lines 
1-2). If the union of all conflict sets is not empty, the search "jumps back" to the 
most recent assignment that is involved in the detected dead end. This means that 
the assignment to the variable Vh is involved in the reached dead end, where h is 
the largest index in this union (cf. Figure |SJ line 3). Before jumping back, the not 
yet defined conflict set of the value assigned to the variable Vh (cf. Figured line 4) 
becomes the union of the conflict sets of the values attempted for the variable V, 
without the index h (Figure |SJ lines 4-5). This is crucial because without any 
change in the assignments to the variables indicated in this conflict set, the variables 
from Vh to Vi will be bound to the same values leading to the same dead end, and 
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thus into a loop. Then the assignments to the variables from Vh to Vi-\ are deleted 
(Figure [3J line 6) and the conflict sets of all previously labelled variables (Vh to V n ) 
are updated, i.e. all defined conflict sets indicating variables that are not deleted 
are kept because they are still valid (Figure |S1 lines 7-15). Finally, h, the index of 
the next variable to be labelled, is returned (Figure line 16). 

The method proposed here is in several respects more general than the orig- 
inal CBJ or its extensions with forward checking (FC-CBJ) also presented 
in l|Prosser 1993|) . or with maintaining arc consistency (MAC-CBJ) presented 
in (Prosser 1995): 

Firstly, our algorithm is not restricted to binary constraints; it processes con- 
straints of arbitrary arities. Secondly, instead of checking each assignment of the 
current variable against the assignments to the already bound variables to deter- 
mine the conflict sets as in the original CBJ, constraint propagation is used in 
our approach to detect inconsistencies and their explanations. This is similar to 
MAC-CBJ (Prosse r 1995|l . where constraint propagation performs arc consistency. 
However, the underlying CHR solver is able to perform stronger, more "global" 
propagation because multi-headed rules allow reasoning over combinations of sev- 
eral constraints: 

Example 4 The single-headed rules of the Boolean CHR solver intro- 
duced in Section perform local propagation maintaining local consistency 
(cf. {Marriott and Stuckey 199^ )), which is the canonical extension of arc consis- 
tency to non-binary constraint problems. Furthermore, the two-headed rules perform 
additional propagation: 

Given: the Boolean variables U, V, X, and Y with domains {0, 1} as well as the 
constraints or (X , U , V) , neg (Y , U) , and or (X , Y , V) . In the first search step, we label 
X = 0. The original CBJ is unable to perform at all because all constraints have un- 
bound variables. Neither forward checking in FC-CBJ nor MAC-CBJ will restrict 
any domains of the not-yet-labelled variables. However, in our approach this la- 
belling triggers the rule or(0,U,V) <=> U=V. This simplifies the constraints to U = 
V. neg(Y,V) and or (X, Y,U) . The equation U=V further triggers the rule neg(Y, V) , 
or(X,Y,V) <=> X=l, Y=0, V=l resulting in an inconsistency. 

Thus, in our approach the assignment to the current variable is not only checked 
against past variable assignments but also against the constraints with future vari- 
ables, maintaining some kind of consistency that is in general stronger than local 
consistency. 

As aforementioned, the conflict sets in our approach are not only stored for 
each variable, they are stored for all possible values of each variable also proposed 
by pSruynooghe 2004| ). This allows us to avoid already detected conflicts after any 
back-jumps or re- assignments to variables: 

Example 5 Let us assume that the value of the variable Vj is in conflict with 
the assignments to the variables V\ and V3 and the value 1 of the variable V? is in 
conflict with the assignments to the variables V2 and V4 . Then, CBJ jumps back to 



14 



Armin Wolf 



int cbj Label (int i) { 

(01) if (var [i] .isBoundO) 

(02) return i+lj 

(03) for (int k=0; k <= 1; k++) { 

(04) if (var[i] . conf lictSet [k] == null) { 

(05) cs . equal (var [i] , k, {i}) ; 

(06) if ( cs , isConsistent () ) 

(07) return i+1; 

(08) else { 

(09) var [i] . conf lictSet [k] = cs .getExplanationO \{i} ; 

(10) cs.delete({i}) ; 

(11) } 

(12) } 

(13) } 

(14) return i; 
} 



Fig. 7. The labelling method of conflict-directed backjumping 



the variable V4, undoing the assignments to the variables V7, Vq, V§ and V4. The 
value recently assigned to the variable V4 is thus known to be in conflict with the 
variables V\, V2 and V3. If there is any non- conflicting assignment to the variable 
V4, we know for future labelling that the value of the variable Vj is still in conflict 
with the assignments to the variables Vi and V3. □ 

A generalisation of this search process for arbitrary finite domains is quite simple: 
It must be iterated over the values in the domain (Figure0 lines 3-13 and Figure|Hl 
lines 9-15), and the tests and calculations must be done for all conflict sets of the 
domain values ( Figure |H1 lines 1 and 3-5). 




5.5 Dynamic Backtracking 

Dynamic backtracking (DBT) ( |Ginsberg 1993| ) is 
a guided depth-first search dynamically changing 
the tree structure, i.e. the variable ordering, which 
goes back to the most recent variable assignment 
that is in conflict with the current variable retaining 
the intermediate assignments (cf. Figure |5J) . DBT 
maintains an elimination explanation for each value 
of every variable. Initially, these sets are not de- 
fined, i.e. null. In DBT, two global variable lists 
are maintained to manage the dynamic changes of 
the value ordering. The list unlabelledVars contains the not-yet-labelled variables, 
while the list labelledVars keeps the already labelled variables. 

If the (last-entered) unlabelled variable is not already bound by constraint pro- 
cessing (Figure lines 1-5) the attempt is made to bind it either to the value 
or 1. This variable is labelled with the first value that is possibly not in conflict 
with other already labelled variables (Figure line 7-8). Thus, the value of the 



Fig. 9. The principle of dy- 
namic backtracking 
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int cbj Unlabel (int i) { 

(01) if (var[i] . conf lictSet [0] == && var [i] . conf lictSet [1] == 0) 

(02) return 0; 

(03) int h = max(var [i] . conf lictSet [0] U var [i] . conf lictSet [1] ) 

(04) var[h] . conf lictSet [var [h] . value ()] 

(05) = (var [i] .conf lictSet [0] U var [i] . conf lictSet [1] ) \{h} ; 

(06) cs. delete (.{h, i - 1}) ; 

(07) for (int j=h; j <= n; j++) { 

(08) for (int k=0; k <= 1, k++) { 

(09) if (var [j] .conf lictSet [k] != null 

(10) && var [j] .conf lictSet [k] !=0 

(11) && max(var[j] .conf lictSet [k] ) >= h) { 

(12) var [j] .conf lictSet [k] =null; 



Fig. 8. The unlabelling method of conflict-directed backjumping keeping formerly 
detected and still valid conflict sets 



global counter is chosen as its unique justification and later stored at the variable if 
no inconsistency arises (Figure ITUl lines 8 and 10). This facilitates any subsequent 
deletion of the assignment and all its consequences computed by the underlying 
Boolean constraint solver (cf. Figure fTOl line 16 and Figure ITTI line 15). 

The search continues with the next unlabelled variable if no inconsistency is de- 
tected (FigureEB lines 10-12). Otherwise, the justifications of the already labelled 
variables that are responsible for the detected inconsistency form the elimination 
explanation of the attempted value (Figure ITUl line 15), the assignment is deleted 
(Figure ITUl line 16), and the next value for this variable is attempted (Figure ITUl 
lines 6-19). If each assignment leads to an inconsistency, a dead end is reached, i.e. 
the variable is added to the list of unlabelled variables and i is returned (Figure ITUl 
lines 20-21), which triggers unlabelling. If the elimination explanation of all values 
of the considered variable are empty, the deletion of all variable assignments will not 
resolve the detected inconsistency. Thus, is returned, indicating the inconsistency 
of the given Boolean CSP (Figure ITTI lines 1-2). If the union of all elimination 
explanations is not empty, the search "goes back" to the most recent assignment 
that is involved in the detected dead end. This means that the assignment to the 
variable bt justified by the maximum h in this union (cf. Figure ^TJ lines 5-12) is 
involved in the reached dead end. Before going back, the elimination explanation of 
the value assigned to the variable bt becomes the union of the elimination explana- 
tions of the values attempted for the most recently tried variable in dbtLabel. This 
causes the justification h (Figure El lines 13-14) to be removed. This is crucial as 
in CBJ (see Section IQJl because otherwise bt will be labelled again with the same 
value leading to the same dead end, and thus into a loop. The assignment to the 
variable bt is then deleted (Figure ^TJ line 15) and added to the unlabelled vari- 



(13) 
(14) 
(15) 
(16) 
} 



} 

} 

} 



return h; 
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ables. Additionally, the elimination explanations of all variables are updated, i.e. 
all defined elimination explanations not containing the justification of the deleted 
assignment are kept because they are still valid (Figure ^2 lines 17-33). 
Example 6 Let us assume that the elimination explanation of the value of the 
variable Vj consists of the justifications of the assignments to the variables V% and 
Vg and that the elimination explanation of the value 1 of the variable Vj consists 
of the justifications of the assignments to the variables V3 and V\\. The most re- 
cent assignment (with the largest justification) is assumed to be to the variable V3. 
Thus, DBT goes back to the variable V3, undoing its assignment. The value recently 
assigned to the variable V3 is thus known to be in conflict with the assignments to 
the variables V2, Vg and Vn. Thus, the elimination explanation of this value is the 
union of the justifications of these variables. If there is any non- conflicting assign- 
ment to the variable V3 , we know for future labelling that the value of the variable 
V7 is still in conflict with the assignments to the variables V2 and Vg . Furthermore, 
any other elimination explanation not containing the justification of the removed 
assignment is still valid. □ 

As a result of the non-chronological constraint deletion, a previously bound vari- 
able that is stored in the list of already labelled variables (cf. Figure [Till lines 2-3) 
may happen to be unbound. Such free variables are filtered out and moved to the 
variables that still have to be labelled (cf. Figure ITTI lines 34-37). 

Finally, the number of the variable that has to be labelled next is returned (Fig- 
ure El line 39). 4 

The method proposed here is in several respects more general than the original 
DBT dGinsberg 19930 or its extensions with forward checking (FC-DBT) or even 
with maintaining arc consistency (MAC-DBT) presented in Ij.Tussien et al. 2 000 ) : 

Firstly, our algorithm is not restricted to binary constraints; it processes con- 
straints of arbitrary arities. Secondly, instead of checking each assignment of the 
current variable against the assignments to the already bound variables to deter- 
mine the elimination explanation as in the original DBT, constraint propagation is 
used in our approach to detect inconsistencies and their explanations. This is simi- 
lar to MAC-DBT <p ussi en et al. 2000|l . where constraint propagation performs arc 
consistency. However, the underlying CHR solver is able to perform stronger, more 
"global" propagation because multi-headed rules allow reasoning over combinations 
of several constraints (cf . Example 0J . 

Unlike other solvers for dynamic CSP <|.Tussien et al. 2 000 ) our underlying adap- 
tive CHR constraint solver which was primarily constructed to solve dynamic CSP, 
is adequate to support DBT for dynamic CSP IjVerfaillie and Schiex 1994j) : justi- 
fications are not restricted to variable assignments; any other constraint may be 
justified, too. Thus, the crucial calculation of elimination explanations performed 
by the method getExplanationO returns the identifiers of the constraints involved 
in the detected inconsistency. 

4 Note that the number is not necessarily the index in the array var because the value ordering 
may change dynamically. 
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int dbtLabel (int i) { 



(01) Variable var = unlabelledVars .removeLast () ; 

(02) if (var.isBoundO) { 

(03) labelledVars . add(var) ; 

(04) return i+1; 

(05) } 

(06) for (int k=0; k <= 1; k++) { 

(07) if (var.elimExpl[k] == null) { 

(08) cs . equal (var , k, {cntr}) ; 

(09) if ( cs.isConsistentQ ) { 

(10) var.num = cntr++; 

(11) labelledVars . add (var) ; 

(12) return i+1; 

(13) } 

(14) else { 

(15) var . elimExpl [k] = cs .getExplanationO \{cntr} ; 

(16) cs .delete ({cntr}) ; 



(17) } 

(18) } 

(19) } 

(20) unlabelledVars . add (var) ; 

(21) return i; 
} 

Fig. 10. The labelling method of dynamic backtracking 

A generalisation of this search process for arbitrary finite domains is quite sim- 
ple: It must be iterated over the values in the domain (Figure lines 6-19 and 
Figure ITT1 lines 19-24 and 28-33), and the tests and calculations must be done for 
all eliminations explanation of the domain values (Figure El lines 2 and 6). 

5.6 "Fancy" Backtracking 

The variant of dynamic backtracking presented in (Baker 1994) - we call it "fancy" 
backtracking - is also implemented and compared to the other search strategies. It 
differs from the original dynamic backtracking only in the unlabelling procedure: 
together with the most recent assignment involved in a detected dead end, all as- 
signments that are directly or indirectly determined by this assignment are deleted, 
too. 

Example 7 Let us assume that the most recent assignment involved in a dead 
end is that to the variable Vj. Further, let us assume that its justification is in 
the elimination explanation of the value of the labelled variable V5. Then, the 
assignment of the value 1 to the variable V5 is determined by the assignment to the 
variable V7. Furthermore, any elimination explanation, e.g. that of the value 1 to 
the variable Vg, containing the justification of the assignment to the variable V5 is 
indirectly determined by the assignment to the variable V7. □ 
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int dbtUnlabeKint i) { 

(01) Variable var = unlabelledVars .removeLast () ; 

(02) if (var.elimExpl[0] == && var . elimExpl [1] == 0) { 

(03) return 0; 

(04) } 

(05) Variable bt ; 

(06) int h = max(var . elimExpl [0] U var . elimExpl [1] ) ; 

(07) for (int j=labelledVars . size ()-l ; j >= 0; j— ) { 

(08) bt = labelledVars.get(j) ; 

(09) if (bt.num == h) { 

(10) labelledVars . remove ( j ) ; 

(11) break; 

(12) } 

(13) bt . elimExpl [bt . value ()] 

(14) = (var. elimExpl [0] U var . elimExpl [1] ) \{h} ; 

(15) cs. delete ({/i}) ; 

(16) unlabelledVars . add(bt) ; 

(17) for (j=0; j < unlabelledVars . size () ; { 

(18) var = unlabelledVars . get (j ) ; 

(19) for (int k=0; k <= 1, k++) { 

(20) if (var. elimExpl [k] != null 

(21) && h e var. elimExpl [k] ) { 

(22) var . elimExpl [k] = null; 

(23) } 

(24) } 

(25) } 

(26) for (j=labelledVars.size()-l; j >=0; j— ) { 

(27) var = labelledVars .get (j ) ; 

(28) for (int k=0; k <= 1, k++) { 

(29) if (var. elimExpl [k] ! = null 

(30) && h G var . elimExpl [k] ) { 

(31) var . elimExpl [k] = null; 

(32) } 

(33) } 

(34) if (Ivar.isBoundO) { 

(35) labelledVars . remove ( j ) ; 

(36) unlabelledVars . add (var) ; 

(37) } 

(38) } 

(39) return labelledVars . size 0+1 ; 
} 



Fig. 11. The unlabelling method of dynamic backtracking keeping formerly detected 
and still valid elimination explanations 

Figure shows our implementation of this proposed extension of dynamic back- 
tracking. The additional loop (Figurc lT21 lines 18-33) calculates the set of all justifi- 
cations of the assignments that must be deleted containing at least the justification 
of the most recent assignment involved in the detected dead end (Figure El mic 
16). Then, all these assignments are deleted (Figure IT21 line 34). Additionally, the 
elimination explanations of all variables are updated, i.e. all defined elimination ex- 



Intelligent search based on adaptive CHR 



19 



planations disjoint to the justifications of the deleted assignments are kept because 
they are still valid (Figure E2 lines 23-25 and 45-47). 

6 Performance Comparison 

We have compared the different search procedures presented in Section together 
with a SICStus Prolog implementation of chronological backtracking based on the 
SICStus Prolog compilation of the CHR handler for Boolean constraints presented 
in Section We applied these procedures, then, to all AIM instances with 50 vari- 
ables. 5 For each instance, the required backtracking or backjumping steps together 
with their elapsed runtime are listed in Tables Q and The runtime was measured 
on a Pentium IV PC with 2.8 GHz running Windows XP Professional, Java 1.4.0 
from Sun 6 and SICStus Prolog 7 3.11.0. In Tabled CBJ means conflict-directed 
backjumping introduced by l|Prosser 1993|) as presented in Section l5~4l DBT means 
dynamic backtracking introduced in ( |Ginsberg 1993| ) while FBT is its "fancy" vari- 
ant introduced in ( Baker 1994 ), both presented in Section IB~5l In Table both 
chronological backtracking procedures - the one presented in Section 15.31 and the 
other implemented in SICStus Prolog - obviously require the same number of search 
steps; they differ only in their runtime. 

Based on these results, we compared the number of search steps and runtime 
in graph form. Figure El shows the "qualitative" comparison, and Figure El the 
"quantitative" comparison. In both figures, the last two groups show the sum- 
mations of the steps/the elapsed runtime for the AIM instances with exactly one 
solution (yes-instances) and for the inconsistent instances (no- instances) . The sum- 
mations show that conflict-directed backjumping performs very well, confirming the 
results in IjProsser 1993): in terms of the number of search steps, conflict-directed 
backjumping (CBJ) requires on average two orders of magnitude less than chrono- 
logical backtracking for all instances and is on average more than one order of 
magnitude faster than all other search procedures, even faster than the SICStus 
Prolog implementation. Furthermore, the performance of the Java and SICStus 
Prolog implementations of chronological backtracking are comparable. Looking at 
the required number of search steps, we find that in a few cases dynamic backtrack- 
ing (DBT) requires marginally more steps than chronological backtracking, which 
is at odds with the statements made in IjBaker 1 994). Surprisingly, its extended 
variant (FBT) proposed in (Baker 1994) often requires more search steps than the 
original version of dynamic backtracking. 

7 Discussion 

The conflict-directed backjumping algorithms presented in fProsser 19931 
IProsser 199 5) compute the conflict sets either from total assignments with respect to 

5 Larger instances tended to take too much time (some over 24 hours). 

6 see http://java.sun.org/ 

7 see http://www.sics.se/sicstus/ 
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int fbtUnlabeKint i) { 

(01) Variable var = unlabelledVars .removeLast () ; 

(02) if (var.elimExpl[0] == && var . elimExpl [1] == 0) { 

(03) return 0; 

(04) } 

(05) Variable bt ; 

(06) int h = max(var . elimExpl [0] U var . elimExpl [1] ) ; 

(07) for (int j=labelledVars . size ()-l ; j >= 0; j— ) { 

(08) bt = labelledVars.get(j) ; 

(09) if (bt.num == h) { 

(10) labelledVars.remove(j) ; 

(11) break; 

(12) } 

(13) bt. elimExpl [bt. value ()] 

(14) = (var. elimExpl [0] U var . elimExpl [1] ) \{h} ; 

(15) unlabelledVars . add (bt) ; 

(16) IntegerSet label = {h} ; 

(17) boolean isChanged = true; 

(18) while (isChanged) { 

(19) isChanged = false; 

(20) for (j=0; j < labelledVars . size() ; j++) { 

(21) var = labelledVars .get (j ) ; 

(22) for (int k=0; k <= 1, k++) { 

(23) if (var. elimExpl [k] != null 

(24) && label n var . elimExpl [k] != { 

(25) var . elimExpl [k] = null; 

(26) label = label U {var.num}; 

(27) labelledVars . remove ( j ) ; 

(28) unlabelledVars. add(var) ; 

(29) isChanged = true; 

(30) } 

(31) } 

(32) } 

(33) } 

(34) cs .delete (label) ; 

(35) for (j=labelledVars.size()-l; j >= 0; j — ) { 

(36) var = labelledVars . get (j ) ; 

(37) if ( Ivar.isBoundO) { 

(38) labelledVars. remove(j) ; 

(39) unlabelledVars . add (var) ; 

(40) } 

(41) } 

(42) for (j=0; j < unlabelledVars . size () ; j++) { 

(43) var = unlabelledVars . get (j ) ; 

(44) for (int k=0; k <= 1, k++) { 

(45) if (var. elimExpl [k] != null 

(46) && var. elimExpl [k] n label != 0) { 

(47) var . elimExpl [k] = null; 

(48) } 

(49) } 

(50) return labelledVars. size 0+1; 
} 



Fig. 12. The alternative unlabclling method of a variant of dynamic backtracking 
deleting additional variable assignments also keeping formerly detected and still 
valid elimination explanations 
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Fig. 13. Comparison of different search strategics on the AIM instances with 50 
variables in terms of their number of backjumping/backtracking steps 
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Fig. 14. Comparison of different search strategies on the AIM instances with 50 
variables in terms of their runtime 
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Table 1. Comparison of different search strategies on the AIM instances with 50 
variables (Part 1) 



Intelligent Backtracking 
AIM instance CBJ a DBT 6 FBT 
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28 
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28 
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24 
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30 
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7 
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7 
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7 
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515612 
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Conflict-directed Backjumping as presented in Section 15.41 
b original version of Dynamic Backtracking I Ginsberg 1993 1 
c extended version of Dynamic Backtracking ( Baker 1994 1 

d backjumping/backtracking steps required to find the solution or detect the inconsistency 
e elapsed runtime on a Pentium IV PC with 2.8 GHz running Windows XP Professional 



the violated constraints or by using forward checking or maintenance of arc consis- 
tency, respectively. The calculation of the elimination explanations in IjBaker 19941 
|Ginsberg 1993|IJussien et al. 2f)00j) during practical experiments was rather similar. 
Thus, we assume that inconsistencies are detected rather late, after a lot of super- 
fluous, unsuccessful assignments, i.e. search steps. In our approach, the underlying 
Boolean constraint solver performs local, but also some "global" constraint propa- 
gation (cf. Example 0J. In general, the search spaces are more restricted. Thus, we 
expect inconsistencies to be detected earlier in the search tree, resulting in more 
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Table 2. Comparison of different search strategies on the AIM instances with 50 
variables (Part 2) 
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39657 
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2 


78 


30 
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916 
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1156 
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28 


438 
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15 


265 


329 
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47 
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7 
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236 
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544654 


254650 


427464 


53 no- instances 


9799110 


3781173 


2977018 



backtracking steps required to find the solution or detect the inconsistency 

b elapsed runtime on a Pentium IV PC with 2.8 GHz running Windows XP Professional 



general conflict sets or elimination explanations and also earlier detection of dead 
ends. 

Example 8 Considering the Boolean constraint solver presented in Section^ and 
the constraints 

neg(X,Y), or(X,Y,Z), neg(Z,U) , 

the application of one of the CHR on the first and second constraint will add the 
syntactical equation Z=l. Further addition of the assignment, i.e. equation U=l, will 
result in an inconsistency, i.e. false, by applying one of the CHR to the actualised 
third constraint, i.e. neg(l,l). Thus, in conflict- directed backjumping and in dy- 
namic backtracking, the assignment U=l is excluded during any further search: the 
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conflict set and elimination explanation of value 1 of the variable U will be the empty 
set. Or, if we consider the constraints 

or(X,Y,Z), neg(Z,l) , 

the assignment X=l triggers a CHR that derives the eguation Z=l, resulting in the 
constraint neg(l , 1) and eventually in false. Now, the assignment X=l is excluded 
during any further search, too. □ 

We assume that this kind of "consistency maintenance" , i.e. constraint propaga- 
tion, is - at least partially - responsible for the absence of the bad performance of dy- 
namic backtracking when applied to 3-SAT problems, as reported in IjBaker 1 994). 

8 Conclusion 

During our review of the adaptive CHR system, we have emphasised the potential 
of this system for explanation-based constraint programming (Jussicn 2001). One 
possibility here is the use of explanations in building explanation-guided search 
algorithms. More specifically, we have demonstrated the simplicity of implement- 
ing sophisticated "intelligent" search strategies in conjunction with a CHR-based 
constraint solver within this system. In this context, 

• "simplicity" means that the implementations are quite straightforward, us- 
ing the interface to the underlying adaptive constraint solver in an obvious 
manner 

• "sophisticated" means that early inconsistency detection accomplished by 
constraint propagation within the underlying solver obviously reduces the 
number of search steps 

Conflict-directed backjumping and dynamic backtracking based on CHR thus gain 
a kind of "consistency maintenance" and the poor performance of dynamic back- 
tracking reported in l|Baker 1994|) does not occur. 

An empirical comparison of the implemented search procedures on the AIM in- 
stances showed that the addition of "intelligence" to the search process may reduce 
the number of search steps dramatically. Even the rather simple conflict-directed 
backjumping strategy outperforms on average all the other strategies tested. Fur- 
thermore, we have shown that the runtime of the Java implementations of the 
intelligent search strategies is in most cases better than the implementations of 
chronological backtracking, even better than the implementation in SICStus Pro- 
log. 

One of the main conclusions in a recent paper on building state-of-the-art 
SAT solvers | |Lynce and Marques-Silva 2002| ) "... is that applying non-chronological 
backtracking is most often crucial in solving real-word instances of SAT." - 
Our implemented "intelligent" search procedures belong to this set of non- 
chronological backtracking solvers. A further development of the presented tech- 
niques ( Midler ~2004j) shows that their performance is comparable to those of these 
state-of-the-art SAT solvers. 
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9 Future Work 

Future work will focus on fast implementation techniques, like counter-based 
or lazy implementations, randomised and heuristic variable selection, and as- 
signment strategies that are commonly used in other state-of-the-art SAT 
solvers jLynce and Marques-Silva 2002| ) as well as the implementation of par- 
tial order dynamic backtracking ( |Ginsberg and McAllester 1994| ) or its generali- 
sation IjBliek 19 98). Furthermore, all these extensions and the presented search 
strategies will be compared with some local search algorithms that might also be 
implemented on the basis of our adaptive CHR system (Wolf 2001a). Further future 
research topics are the implementation of the generalisations proposed during our 
presentation of the search strategies and their application to and comparison with 
other finite-domain constraint satisfaction problems like job-shop scheduling. 

Conflict-directed backjumping performs well for Quantified Boolean Logic Satis- 
fiability ( |Giunchiglia et al. 2001| ). The algorithm presented here might therefore be 
adapted and successfully applied to this problem class, which is strongly related to 
SAT. 
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